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Abstract 

We present an initiatory study of quark helicity flip generalized parton distributions 
(GPDs) in rif = 2 lattice QCD, based on clover-improved Wilson fermions for a large 
number of coupling constants and pion masses. Quark helicity flip GPDs yield essential 
information on the transverse spin structure of the nucleon. In this work, we show first 
results on their lowest moments and dipole masses and study the corresponding chiral 
and continuum extrapolations. 



1 Introduction 

Generalized parton distributions (GPDs) [1] have opened new ways of studying the complex 
interplay of longitudinal momentum and transverse coordinate space [2,3], as well as spin 
and orbital angular momentum degrees of freedom in the nucleon [4]. As a counting of the 
helicity amplitudes in Fig. [T] reveals [5], there are eight independent real functions needed at 
twist 2. Four of them, namely Ht, Et, Ht and Et, are related to a flip of the quark helicity, 
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Figure 1: The lower part of the handbag diagram. 



/i = — //, hence quark helicity flip GPDs 1 . Quark helicity flip GPDs play a prominent 
role in the understanding of the transverse spin structure of the nucleon and significantly 
sharpen positivity bounds on GPDs in impact parameter space [6]. Specifically, it could be 
very interesting to exploit and study the equation-of-motion relations between the lowest 
moments of quark helicity flip, unpolarized and twist-3 GPDs which have been obtained 
in [6]. The (chirally odd) tensor GPDs also provide a framework with which to study the 
correlation between quark spin and quark angular momentum in unpolarized nucleons [7]. 

Quark helicity flip GPDs are defined via the parameterization of an off-forward nucleon 
matrix element of a quark operator involving the a^-tensor as follows [5] 



Here the momentum transfer is given by A = P' — P with t = A 2 , P = (P' + -P)/2, and 
£ = — n - A/2 denotes the longitudinal momentum transfer, where n is a light-like vector. 
The first of these tensor GPDs, H T (x,^,t), is called generalized transversity, because it 
reproduces the transversity distribution in the forward limit, H T (x, 0,0) = 5q(x) = h\{x). 
Integrating HT(x,£,t) over x gives the tensor form factor: 



Since the quark tensor GPDs require a helicity flip of the quarks, they do not contribute 
to the deeply virtual Compton scattering (DVCS) process 7*p — > jp'. Naively, one could 
think that this could be balanced by the production of a transversely polarized vector meson 
instead of a photon, j*p — ■> m T p'. However, it has been shown that the corresponding 
amplitude, remarkably, vanishes at leading twist to all orders in perturbation theory [8-10]. 

1 Also called tensor GPDs. 




E T (x,t,t) 



(1) 
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The only process giving access to the generalized transversity which has been proposed in the 
literature so far is the diffract ive double meson production j*p — > m L m T p' [11]. Naturally, 
one expects the measurement of this reaction to be much more involved than e.g. the exclusive 
electroproduction of a single vector meson. Since the tensor GPDs are practically unknown, 
it is unclear how to even estimate the corresponding cross section to see if a measurement 
of this process is at all feasible. Given that the situation seems to be much more difficult 
than for the (un-)polarized GPDs, lattice calculations of the lowest moments of the quark 
helicity flip GPDs will be highly valuable. While (un-)polarized GPDs have already been 
investigated in a number of papers [12-19], we present here the first lattice calculation of 
quark helicity flip GPDs. 

Lattice calculations of moments of parton distributions mostly disregard the computation- 
ally expensive quark-line disconnected contributions. They correspond to a situation where 
the operator is inserted into a closed quark loop which is connected to the nucleon only via 
gluons. Since the tensor operators flip the quark helicity, these disconnected diagrams do 
not contribute in the continuum theory for vanishing quark masses. Therefore, we expect 
only small contributions for the disconnected graphs in our calculation. This expectation is 
supported by numerical results from [20], where the tensor charge was calculated in quenched 
lattice QCD. The authors explicitly computed the disconnected pieces for the tensor operator 
and found the contributions from up- and down-quarks to be compatible with zero within 
one standard deviation. Thus it is possible to estimate the individual up and down quark 
tensor GPDs, which is a major advantage compared to other observables where usually only 
the iso-vector channel is considered. Further early results on the tensor charge in quenched 
lattice QCD have been presented in [21,22]. 

As mentioned above, in calculating the lowest moments of the tensor GPD H T (x,£,t), we 
automatically obtain the corresponding moments of the transversity distribution, (x™ -1 )^, 
for t = £ = 0. The quark transversity has recently attracted renewed attention related to 
the Collins asymmetry in e.g. semi-inclusive deep inelastic scattering. It is generally believed 
that transverse single-spin asymmetries (SSA) [23] are generated predominantly by the Sivers 
and Collins mechanism. These two differ in their dependence on the azimuthal angles and 
thus can be separated. The contribution due to the Collins mechanism is proportional to 
a convolution of the transversity distribution 8q(x) and the Collins fragmentation function 
Hi(z), which are both chiral odd. Lack of knowledge of both the transversity and the 
Collins function, however, seriously hampers the interpretation of the exciting experimental 
results on such SSAs [24,25]. Lattice results for the lowest moments of 5q(x) for up and 
down quarks could help to reveal the physics behind these measured asymmetries. 

The paper is organized as follows. We begin by briefly reminding the reader of the methods 
and techniques we use to extract moments of GPDs from the lattice in Section 2. In Section 3, 
we specify the parameters of our calculation and present our results for the lowest moments 
of the tensor GPD Ht(x, £, t). Making use of the large number of results for different sets of 
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lattice parameters, we attempt to extrapolate the moments of the generalized transversity 
as well as the dipole masses of the tensor GPDs to the continuum and chiral limits. Finally, 
in Section 4 we summarize our findings. 



2 Extracting moments of GPDs from lattice simula- 
tions 

On the lattice, it is not possible to deal directly with matrix elements of bi-local light-cone 
operators. Therefore, we first transform the LHS of Eq. to Mellin space by integrating 
over x, i.e. J_ 1 dxx n ~ 1 . This results in nucleon matrix elements of towers of local tensor 
operators 

O^ 1 -^- 1 (0) = ?p{0)tcT^ u tD^ . . . W^- l} iP{0) , (3) 

which are in turn parameterized in terms of the tensor generalized form factors (GFFs) 

Axni, Bxni, ^-Tni and Bmi- Here and in the following, D = \{D — D) and {• ■ • } indicates 
symmetrization of indices and subtraction of traces. The parameterization for arbitrary n is 
given in [26,27] 2 . Here we show explicitly only the expressions for the lowest two moments. 
For n = lwe have 

(P'A'\ ^(0) W "^(0) |PA) = U{P\ A') {ia^A T10 (t) + *—^-A T10 (t) 

+ 1 ^B TlQ {t)\U{P,K) . (4) 

The inclusion of an additional term oc j\pp"^ = ^p v — ^ V P^ in Eq. is forbidden by time 
reversal symmetry [5]. For n = 2, however, this can be balanced by including another factor 
of A, leading to four generalized form factors, 



A [H S {u , l} (P'A'I ^(0)^1^(0) \PA) = A [H S { ^ l} U(P', A') {ia** 'P w A T20 (t) 

+ ^~P^A T2Q {t) + y^-P^B^t) 



+ A^B T21 (t)}U(P,A), (5) 

m 

up to trace terms, where A^i and denote anti-symmetrization and symmetrization 

of (n, z/), respectively. For n = 3 there are seven independent tensor GFFs, as an explicit 

2 Note that the Mellin-moment index n used here differs from the number of covariant derivatives n in [26] 
by one. 
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counting shows [26,27]. The simultaneous extraction of such a large number of GFFs poses 
a challenge for lattice QCD calculations, which we plan to address in the near future. 

Instead of calculating continuum Minkowski space-time matrix elements (e.g. in Eqs. (J3J) 
and (|HJ)) directly, on the lattice we work within a discretized Euclidean space-time framework 
to calculate nucleon two- and three-point correlation functions. The nucleon two- and three- 
point functions are given by 

C 2pt (r,P) = ^TikiN^P^r^P)), 
C 3 ^-^(t,P',P) = ^f ifc (iV fc (T s ^,PO0r W ^"" 1 (^i(^c,P)>, (6) 

j,k 

where T is a (spin-)projection matrix and the operators N and N create and destroy states 
with the quantum numbers of the nucleon, respectively. The relation of C^f* to the parame- 
terizations in Eqs. (jU and (jSJ) is seen by rewriting Eq. (JHJ) using complete sets of states and 
the time evolution operator, 

r ,3ptAt I ^i...Mn-l/_ p> p\ _ { Z (P) Z (P )) E(P)( T -T stc )-E{P')(T mk -T) 

lT '^'^" 4P(P)P(P) 

x ( p '> A 'l Or" l-P, A) U(P, A)VU(P', A') + . . . . (7) 

A, A' 

Similarly, the two-point function for T = 1/2(1 + 74) can be written as 

C 2 P\T, P) = (Z(P)Z(P)) 1/2 E ^ + m e -E{P){r-r^) + _ _ (g) 

P(P) 

The ellipsis in Eq. (JZJ and (jHJ) represents excited states with energies E' > E(P), E(P'), 
which are exponentially suppressed as long as r — r src ^> l/P",T sn k — r 3> 1/E'. Inserting 
the explicit parameterizations from Eqs. (jlj) and (jHJ transformed to Euclidean space into 
Eq. (jjj), we sum over polarizations to obtain 

(Z(P)Z(P')) 7 

. -i< i')lr-7\,, l-Ai / ' )( r„ 1 .. I p^p' _ rfl j 



^Spt/i^i-Mn-l/ p/ pN = V^V-t -j?(P)(r-r 5rc )-i;(P / )(r snk -r) T 

K ' ' y AE(P')E(P) 



x C-^^W + ^-^B Tn0 {t) + •••){#> -m) 



(9) 



where e.g. a^ 41 is the Euclidean version of the prefactor ia^P^ 1 in Eq. (|5|). The Dirac- 
trace in Eq. is evaluated explicitly, while the normalization factor and the exponentials 
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in Eq. (JJJ) are cancelled out by constructing an appropriate r-independent ratio R of two- 
and three-point functions, 



Ro(r,P',P) 



^ pt (r,P',P) r c^jr, P')C^(r snk , P')C^(r SQk - r + r sl . c , P) l * 
C^H^nk, P') [ C**(t, P)C^(r snk , P)C^(r snk - r + r src , P') _ 



(10) 



The ratio R is evaluated numerically and then equated with the corresponding sum of 
GFFs times P- and P'-dependent calculable pre-factors, coming from the traces in Eq. Q. 
For a given moment n, this is done simultaneously for all contributing index combinations 
(fjLVfjLi . . . Hn-i) and all discrete lattice momenta P, P' corresponding to the same value of 
t = (P' — P) 2 . This procedure leads, in general, to an overdetermined set of equations from 
which we finally extract the GFFs [16]. We have taken care to ensure that our normaliza- 
tion leads exactly to the x-moment of the transversity distribution 5q(x) = h\{x) as defined 
in [28] . To make this as transparent as possible, we give an explicit example of one of the 
equations we use to extract (x)s 



where only the Ti (see Eq. fJ17jl ) projector contributes and 2{34} represents the operator 



On the lattice the space-time symmetry is reduced to the hypercubic group H(4), and the 
lattice operators have to be chosen such that they belong to irreducible multiplets under 
H(4). Furthermore, one would like to avoid mixing under renormalization as far as possible. 
In the case of the twist-2 operators in Eq. or more precisely their Euclidean counterparts, 
this presents no problem for n — 1 and n — 2, the only cases to be considered in this paper. 
For n = 1 we have the 6-dimensional multiplet consisting of the operators 



which is irreducible in the continuum as well as on the lattice {H(A) representation r{ in 
the notation of [29]). The 16- dimensional space of continuum twist-2 operators with n = 2 
decomposes into two 8-dimensional multiplets transforming according to the inequivalent 
representations t[ 8 ^ and r!f\ Typical members of these multiplets are, e.g., 




(11) 



^(0)i<vV(0) , 



(12) 




(13) 



in the case of r{ , and 




(14) 
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Table 1: Lattice parameters: gauge coupling j3, sea quark hopping parameter K sea , lattice 
volume, number of trajectories, lattice spacing and pion mass. 

(s) 

for r 2 . All these operators are free of mixing problems, but one has to take into account that 
operators belonging to inequivalent representations have different renormalization factors. 

Obviously, for a successful computation of the GFFs, one would like to have as many 
different nucleon sink and source momenta and projection operators T as possible in order 
to obtain a large number of independent non- vanishing Dirac-traces in Eq. (0). This is 
particularly true for the tensor operators because they involve a^ v and the number of tensor 
GFFs grows rapidly with n. Once we have extracted the GFFs from the lattice correlation 
functions, it is an easy exercise to reconstruct the corresponding moments of tensor GPDs, 
H^(^,t) = f dxx n ~ 1 H T (x,£,,t) etc., using the polynomiality relations [26] 

Hr\U)=A TW (t)=g T (t), 
Hr\U)=A TW (t), 

Er 1 (^t) = B T10 (t), 

Er 1 (Z,t) = B T10 (t) = 0, 

These equations directly show that for n < 2, a dependence on the longitudinal momentum 
transfer £ is only seen for the GPD Et, which is the only quark GPD odd in £. In order 
to investigate the £ dependence of the generalized transversity H^(^,t), one has to consider 
at least the n = 3 Mellin moment. Finally, we note that in the forward limit the moments 
H^i^t) reduce to the moments of the transversity distribution, if^(£ = 0,t = 0) = (x™ -1 )^ 



Tjn=2 



rjn- 

rirp 
rpn= 

Jllrp 

rpn=2 

Ulrp 



(Z,t)=A T20 (t), 
■ 2 (U) = A T20 (t), 
2 (£,t) = B TW (t), 
(£,*) = (-20Srai(t). 



(15) 



7 




Figure 2: The tensor form factor together with a dipole fit. 




-t[GeV 2 ] -t[GeV 2 ] 
Figure 3: The GFF At2q together with a dipole fit. 



3 Lattice results for moments of the generalized 
transversity 

The simulations are done with nf = 2 flavors of dynamical non-perturbatively 0(a) improved 
Wilson fermions and Wilson glue. For four different values (3 = 5.20, 5.25, 5.29, 5.40 and 
three different k values per j3 we have in collaboration with UKQCD generated (9(2000 — 
8000) trajectories. Lattice spacings and spatial volumes vary between 0.07-0.11 fm and 
(1.4-2.0 fm) 3 respectively. A summary of the parameter space spanned by our dynamical 
configurations can be found in Table ^ We set the scale via the force parameter r , with 
r = 0.467 fm. 

Correlation functions are calculated on configurations taken at a distance of 5-10 trajecto- 
ries using 4-8 different locations of the fermion source. We use binning to obtain an effective 
distance of 20 trajectories. The size of the bins has little effect on the error, which indicates 
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Figure 4: The effective dipole mass as a function of a cut in t. 
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auto-correlations are small. In this work, we simulate with three choices of sink momenta 
P' and polarization operators, namely 



2n - 2ir 

P^ = (0, 0, 0), Pi = ( — , 0, 0), P' 2 = (0, — , 0) , 



where Ls is the spatial extent of the lattice, and 

1 unpol 



-(1 + 74), f 1 = ^(1 + 74) «*757l> r 2 = -Al + 74) il5l2 ■ 



(16) 



(17) 



The choice of the two polarization projectors, f 1 and T 2 is particularly advantageous for the 
extraction of the tensor GFFs. The values of the momentum transfer A = (2n/Ls) q used 
in this analysis are 

q: (0,0,0), (1,0,0), (1,1,0) (1,1,1), (2,0,0) (18) 

and the vectors with permuted components. All lattice results below have been non- 
perturbatively renormalized [30] and transformed to the MS scheme at a renormalization 
scale of 4 GeV 2 . 

In this work, we focus on the lowest two moments of the GPD Ht- A broader analysis will 
in particular include moments of the linear combination 2Ht{x, £, t) +Et{x, £, t) which have 
been shown to play a fundamental role for the transverse spin structure of the nucleon [6]. 
Furthermore, in [7] it is claimed that the x- moment of this linear combination gives the 
angular momentum carried by quarks with transverse spin in an unpolarized nucleon, in 
analogy to Ji's sum rule. In Figs. El and E2 we show our results for the lowest two moments 
of the generalized transversity for up and down quarks in the nucleon as functions of the 
squared momentum transfer t. The lattice points and dipole curves are the result of a 
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combined dipole fit together with linear continuum and pion-mass extrapolations of the 
form 



4 dipole,m w ,a/ + \ _ ^Tno(O) + a l m l + a 2Q- 2 



with five fit parameters ^4^(0), m° D and ai,...,as. The curves show the fit function 
in the continuum limit, i.e. for a = 0, at the physical pion mass. Correspondingly, the 
difference A^^ e)m " K ' a (t) — j\^P° le ' m ^ ,a_0 (t) has been subtracted from the individual data 
points before plotting. Although the extrapolation to the continuum limit turns out to be 
almost flat, except for A^O) for which a 2 ~ —4.2 ± .7 fm~ 2 , we include the a 2 -dependence 
because it reduces the x 2 of the fits considerably. To check our ansatz in Eq. (jl9j) . we show in 
Fig-IUthe (effective) dipole mass m° D as a function of a cut for minimal and maximal values of 
the momentum transfer squared t used for the fit, t min < t < t max (keep in mind that t < 0). 
The effective dipole mass is in both cases very stable and constant, except when — t max 
becomes large since there are not enough data points used in the fit to determine the dipole 
mass accurately. Still, a more sophisticated approach is desired for future investigations. 
Additionally, the assumed linear dependence on a 2 and m 2 eventually has to be replaced by a 
functional form obtained from e.g. chiral perturbation theory. The quark mass dependence 
of the first two moments of the (iso-vector) transversity has already been investigated in 
[31,32]. 

The forward moments and dipole masses at m n = mf ys and a = are found to be 

= A% 10 (0) = 0.857±.013, m D = 1.732±.036 GeV, 

(1)1 = A% 10 (0) = -0.212±.005, m D = 1.741±.056 GeV, 

(x)Jf = A" 20 (0) = 0.268±.006, m D = 2.312±.071 GeV, (20) 

(x) d s = Af 20 (0) = -0.052±.002, m D = 2.448±.173 GeV, 



and for the iso-vector and iso-singlet combinations we obtain the dipole masses 
A T io ■ m u D d = 1.731±.034 GeV, m u + d = 1.713±.043 GeV, 

A T 2o ■ m u ~ d = 2.318±.067 GeV, m u + d = 2.286±.083 GeV, ^ ' 

which agree with the up- and down-quark dipole masses within errors. Our result for the 
iso-vector tensor charge (l)^~ d = 1.068 ± 0.016 is in agreement with results in [20] and 5% 
to 15% lower compared to lattice studies in [32-35]. However our result for the iso-vector 
x-moment (x)^~ d = 0.322 ± 0.006 is substantially lower than the quoted value of (x)$~ d = 
0.533 ± 0.083 (unquenched, n = 0.1570, from [34]) and also the chirally extrapolated value 
(x) u f d = 0.506 ± 0.089 [32] 3 . Since previous works used unimproved Wilson fermions with 

3 This holds also for up and down quarks separately. 
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Figure 5: The pion mass dependence of the lowest two moments of the transversity distri- 
bution. 



no continuum extrapolation together with perturbative renormalization of the operators, 
the numbers should be compared with some care. Still, the discrepancy could indicate some 
problems with the normalization. 

The explicit dependence of the tensor charge gxit = 0) = (1)^ and the x-moment of the 
transversity (x)$ on the pion mass is shown in Fig. where all points have already been 
extrapolated to the continuum limit. The linearly extrapolated values at mf ,s agree within 
errors with the results from the global fit in Eq. (f2T)|) . From the figures we see that the tensor 
charge is approximately constant over the available range of pion masses, while e.g. (x)$ 
clearly shows a dependence on m^phys and drops by ~ 20% going from m 2 = 1.4 GeV 2 down 
to ml = 0.4 GeV 2 . 

Interestingly, our results for the iso-vector dipole masses for the first two moments of 
agree very well with those obtained from fits to the moments of the polarized GPD, H [36], 
which are shown to lie on a linear Regge trajectory. It will be interesting to see if this trend 
continues for higher moments. 
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Figure 6: The ratio in Eq. (|2*3|) as a function of m^. for n — 0, 1. 



Finally, in Fig. El we investigate the Soffer bound [37] 



<-(Aq(x)+q(x)) 



(22) 



which holds exactly only for quark and anti-quark distributions separately. Mellin moments 
of the distribution functions as defined in section |21 give however always sums/differences of 
moments of quark and anti-quark distributions, e.g. (x n ) q + (— l) n+1 (x n )q. Taking Mellin 
moments of Eq. (}2*2*|) and assuming that the antiquark contributions are small, we expect 
that the ratio 



2\(x n )i 



((x n ) + (x n ) A ) 



n 



0,1 



(23) 



is smaller than one. In Fig. El we show this ratio for up and down contributions as a function 
of m^. As we can clearly see from the figure, the ratio in Eq. ()23|) is smaller than one over 
the whole range of available pion masses. Taking into account what has been said above, this 
strongly indicates that the Soffer bound is satisfied in our lattice calculation of the lowest 
two moments of the unpolarized, polarized and transversity quark distributions. 
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4 Conclusions and Outlook 



We have computed the lowest moments of the quark tensor GPD Ht in lattice QCD and 
studied the chiral and continuum limit of the forward moments and the dipole masses. 
Assuming that contributions from anti-quarks are small, our results indicate that the Soffer 
bound, relating the transversity, unpolarized and polarized quark distributions, is satisfied 
in our calculation. 

The results are promising and our study will soon be extended to include the tensor GPDs 
Et, Ht and Et- Once a set of the lowest moments of all tensor GPDs is available, it will 
be extremely interesting to analyze the transverse spin density of quarks in the nucleon, 
the corresponding positivity bounds and the relation to moments of twist-3 GPDs using 
sum-rules obtained from the equation of motion [6] . 
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